Phason Disorder Effects in the Penrose Tiling Antiferromagnet 
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We discuss the ground state of a disordered two dimensional Heisenberg antiferromagnet. The 
starting structure is taken to be a perfectly deterministic quasiperiodic tiling, and the type of disorder 
we consider is geometric, involving frozen phason flips of a randomly selected subset of sites. We 
consider S=l/2 quantum spins placed on the vertices of the tiling, and interacting with the nearest 
neighbor spins with a uniform exchange interaction J. We calculate the energy spectrum, ground 
state energy and real space local magnetization values as a function of degree of disorder. We find 
that quantum fluctuations are enhanced by disorder. The real space staggered magnetization loses 
its symmetry properties and the average staggered magnetization decreases compared to the case 
of the perfect Penrose tiling. We explain our results in terms of a simple Heisenberg star cluster 
model that takes into account the changes of local environments due to phason flips. 

PACS numbers: 71.23.Ft, 75.10.Jm, 75.10.-b 



It is "well-known that the electronic properties of ran- 
dom tilings are qualitatively different from those of the 
perfect quasiperiodic tiling To cite one example, the 
perfect quasiperiodic: system tends to have strongly sin- 
gular fluctuations in the density of states, while the ran- 
dom systems have fluctuations that are similar in some 
respects to those present in weakly disordered periodic 
systems. We consider a disordered magnetic model : 
the Heisenberg antiferromagnet defined in disordered ap- 
proximants of the two dimensional Penrose tiling 0, Q • 
We investigate the excitation spectrum, the ground state 
energy and staggered local magnetizations in these sys- 
tems and compare the results with the perfect determin- 
istic case (see We create a frozen disorder in the 
perfect deterministic tiling by making randomly located 
phason flips. Phason flips are operations that reorga- 
nize the structure locally, in the vicinity of the flipped 
site. This type of disorder is strongly constrained, as one 
does not modify the basic building blocks of the struc- 
ture, but only the way they are put together. Many 
structurally refined quasicrystals are assumed to be ba- 
sically deterministic, with phonon and phason disorder 
at finite temperatures [2|. However, models for basi- 
cally random structures are also proposed for a second 
category of quasicrystals, based on arguments related to 
maximising entropy of the system Q . Many experimen- 
tal studies, using neutron or X-ray diffraction have been 
carried out to detect phason disorder and compare with 
existing theoretical models. There are many studies of 
physical properties of random quasiperiodic structures 
and the influence of phason disorder. For example, it 
is interesting to consider the effect of disorder on trans- 
port properties, where some studies have shown that the 
electrical conductivity improves when disorder increases. 



The Hamiltonian We consider the antiferromagnetic 
spin-i Heisenberg model with nearest-neighbor interac- 



tions, 



H 



Sj, J>0, 



(1) 



where are pairs of sites that are linked by an edge 
of the tiling. Since the structures we consider are con- 
structed from four sided polygons (the thick and the thin 
rhombus), the model is bipartite (i.e. can be divided 
into two sublattices, A and B). The antiferromagnetic 
couplings are thus unfrustrated. At T=0, since there is 
no frustration in the model, one expects a Neel ordered 
ground state, with equal and oppositely directed sub- 
lattice magnetizations. This ground state expected to 
slightly deviate from the ground state of classical spins 
on the same vertices. The ground state of classical spins 
is (Si) = (+) on a sublattice and (-) on the another 
sublattice. This deviation is expected to be small and 
one could calculate this deviation using linear spin wave 
theory (@, [§] and references therein). At finite temper- 
ature, the long range order will be destroyed due to the 
Mermin-Wagner-Hohenberg theorem , although short 
range correlations will persist. In our calculations we con- 
sider rectangular (Taylor) approximants of the Penrose 
tiling. These are finite samples of N spins satisfying peri- 
odic boundary conditions. Approximants corresponding 
to the perfect deterministic Penrose structure were first 
obtained by a standard method and then these were ran- 
domized using the procedure described below. 

Phason-flip A "phason flip" is a process in which a 
3-fold site hops to a new allowed position. In the process, 
the three bonds linking it to its original nearest neigh- 
bors are effaced, while new connections appear to three 
new nearest neighbors (Fig. [T|). Note that the original 
site and the new phason shifted site belong to different 
sublattices. To maintain a ground state of total spin 
zero, our random selection of phason flip centers is done 
with the requirement that there be an equal number of 
flips on each of the sublattices A and B. In our random 
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FIG. 1: A single phason flip showing the original (ro) and 
final positions (r' ). The original sites and bonds are shown 
in blue and the new site and new bonds in red. 



tiling generating procedure, each step randomly selects 
one three-fold site from a list of possible candidates for 
phason flips. The new connectivity matrix is determined 
and the set of sites available for the next phason flip is 
determined. For measuring the disorder in a given value 
of phason flips, we calculating the overlap between the 
original and the phason disordered tiling. We calculate 
the number of shifted sites and this number normalized 
with the system size will be the degree of disorder (A) for 
a given number of phason flips. The system sizes that we 
considered are N=644, 1686 and 4414 spins. The number 
of realizations considered range from 100 for the smallest 
size, to 10 realizations in the biggest size. 
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have made a study of the evolution of the frequency of 
sites of each coordination number z — 3, 10 as a func- 
tion of the degree of disorder. As shown in Fig. [2l the 
disorder clearly increase the number of 4-fold sites at the 
expense of the 3- fold and the 5-fold sites. The original 
perfect structure has a large number of 3-fold and 5-fold 
sites, with the most common cluster found in the origi- 
nal perfect Penrose tiling, the "football" -shape five- fold 
stars. These clusters are destroyed progressively by dis- 
order. The evolution of coordination numbers will be sig- 
nificant for the qualitative analysis that we will present 
below for the staggered magnetization. 

Linear spin wave analysis Once the randomized 
samples are obtained, a linear spin wave analysis is car- 
ried out as described in [H H|. The A and B sublat- 
tice spin operators are transformed using the Holstein- 
Primakoff transformation to boson operators aj , bj — 
1, N/2). The linearized Hamiltonian in the boson op- 
erators was then diagonalized numerically. Once the 
eigenmodes and spectrum have been determined, one can 
find the ground state energy, and local staggered magne- 
tizations for each realization of disorder, and carry out 
the statistical analyses of the results. We now describe 
the results for the ground state energy, the energy spec- 
trum and local staggered magnetizations. 
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FIG. 3: Density of states plotted in the case of different degree 
of disorder: A (The system size is N=4414). 



FIG. 2: Dependence of the coordination number distribution 
from the degree of disorder (A). On the y axis: the ratio of 
the number of the sites with a given coordination number per 
the system size in the case of certain degree of disorder. (The 
system size we considered is N=4414.) 

Geometrical effects Phason flips clearly lead to a 
modification of local environments on the tiling, and in 
particular, the apparition of new local environments that 
were not present in the original tiling such as 8-fold, 9- 
fold and 10-fold vertices. Furthermore, finite-sized re- 
gions of perfect five-fold rotational symmetry progres- 
sively disappear with increasing degree of disorder. We 



Energy spectrum We show the density of states 
(DOS) curves as a function of the energy is expressed in 
the unit of J for systems with different degree of disorders 
(Fig. [3]). The primary effect of increased phason disorder 
is to smoothen the DOS curve and fill in the gaps. The 
low energy tail of the integrated curve is quadratic, and 
can be fitted to give a spin wave velocity that increases 
with the degree of disorder (Fig. |4j). This indicates that 
spin wave propagation is facilitated by the phason dis- 
order, in analogy with the problem of quantum diffusion 
of electrons in the tight binding model in quasiperiodic 
tilings [1]. Changes in the DOS curve are quite small 
below the degenerate peak E = 3 - disorder only weakly 
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FIG. 4: Spin wave velocity as a function of disorder (A). The 
system sizes we considered are 644, 1686 and 4414. 
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FIG. 6: Comparison of the real space local staggered magneti- 
zation values on the Penrose tiling approximant. We colorized 
the sites on the Penrose tiling using different colors for dif- 
ferent local staggered magnetization value. The color code is 
the following: red, orange, yellow, green, blue with local stag- 
gered magnetization from the lowest to the highest values. 
The system size is 1686. 



affect the "long wavelength" wavefunctions. The degen- 
erate peak at E = 3 is reduced due to the progressive 
disappearance of 3-fold sites, as mentioned already. The 
DOS increases with degree of disorder between energies 
3.4 and 4.3, which is a range of energy corresponding to 
wavefunctions located mostly on four-fold sites. 
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FIG. 5: Ground state energy as a function of increasing dis- 
order (A). The straight line is the ground state energy for 
the perfect Penrose tiling. 



Ground state energy Fig. [5] shows the decrease of 
the ground state energy per site Eq as a function of in- 
creasing degree of disorder. The decrease of Eq is similar 
to the decrease in the average value of the staggered mag- 
netization (see below), showing that the quantum fluctu- 
ations are increased in the random tiling as compared to 
the perfect case. 

Local staggered magnetizations The real space 
staggered magnetizations for a randomized system m S i — 
\{Si)\ are shown in a color plot on the right hand side 
figure in (Fig.[6|). The right hand figure shows the corre- 
sponding quantities on the perfect tiling, for comparison. 
The perfect case shows symmetries that the random case 
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FIG. 7: Average local staggered magnetization in the case 
of different degree of disorder (A). The straight line is the 
average local staggered magnetization for the perfect Penrose 
tiling approximant. We considered the expectation values for 
the thermodynamic limit. In the inset: local staggered mag- 
netization on the central site of a two level Heisenberg-star 
cluster. We use the same coordination number distribution 
like the whole tiling has in the case of different degree of dis- 
order (A). The straight line is the average local staggered 
magnetization for the two level Heisenberg-star cluster using 
the same coordination number distribution like the perfect 
Penrose tiling. 



does not. In addition, the overall (site-averaged) stag- 
gered magnetization is clearly smaller for the disordered 
case (the color coding shows a redder tone, indicating 
smaller values). Fig. [7] shows the average staggered mag- 
netization as a function of degree of disorder. This quan- 
tity decreases with increasing disorder, indicating greater 
quantum fluctuations. 

The results for the staggered magnetizations can be 
qualitatively explained in terms of a finite cluster model 
[5|. In this approximation, the local environment is re- 
tained up to second neighbor level in the spin wave cal- 
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culation. Boundary conditions, and connectivity beyond 
second neighbors are thus neglected. Nevertheless, with 
this simple approximation many of the results can be ob- 
tained, using data about nearest neighbor and next near- 
est neighbor coordination number as input structural pa- 
rameters. We conclude from this that the local staggered 
magnetization values are dependent mainly on the local 
environments for the perfect as well as the disordered 
Penrose tilings (Fig. [7] inset). 

Conclusion We investigated the effects of randomly 
located "phason" flips on a Heisenberg antiferromagnet 
in approximants of the Penrose tiling. We calculated the 
energy spectrum, the ground state energy and local stag- 
gered magnetizations in the ground state using linear spin 
wave approximation in order to investigate the effects of 
increased phason disorder. The magnon speed estimated 
from the low energy part of the spectrum increases with 
phason disorder. The ground state energy is observed to 
decrease with increasing disorder. We investigated the 
changes of the staggered magnetizations in real space. 
The ground state loses its symmetry properties and the 
overall staggered magnetization decrease with increasing 
degree of disorder. This observation, along with the de- 
crease of the ground state energy, show that disorder 
tends to increase quantum fluctuations in this system. 
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